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Abstract. - Drops deposited on rough and hydrophobic surfaces can stay suspended with gas 
pockets underneath the liquid, then showing very low hydrodynamic resistance. When this su- 
perhydrophobic state breaks down, the subsequent wetting process can show different dynamical 
properties. A suitable choice of the geometry can make the wetting front propagate in a stepwise 
manner leading to square-shaped wetted area: the front propagation is slow and the patterned sur- 
face fills by rows through a zipping mechanism. The multiple time scale scenario of this wetting 
process is experimentally characterized and compared to numerical simulations. 



When a liquid droplet is deposited on hydrophobic 
micro-structured materials, it can bead off and stay sus- 
pended with very high contact angles (Cassie-Baxter state, 
hereafter CB) and very low hydrodynamic resistance ("Lo- 
tus effect" ) [1-4] . This property makes such materials use- 
ful for a wide series of applications ranging from coating to 
microfluidics [5-7]. However, when superhydrophobicity 
breaks down, fluid can enter and fill the micro-structures, 
resulting into a smaller effective contact angle (Wenzel 
state, hereafter W). In some situations this transition from 
the CB to the W state is highly desirable. An example is 
provided by heterogeneous porous catalysts, where super- 
hydrophobicity is an unwanted effect as it reduces the con- 
tact area [8]. Understanding the mechanisms triggering 
the transition and characterizing its dynamical properties 
is crucial. The triggering mechanisms for the transitions 
are still widely debated in the literature [9] . This is due to 
the presence of energy barriers encountered when passing 
from CB to W that may strongly depend on the filling pro- 
cedure. This means that even when the Wenzel state has 
lower energy (lower contact angle), the transition is not al- 
ways spontaneous and has to be triggered, for example, by 
pressing on the drop [3] or by using volumetric forces such 
as gravity [9]. However, under certain conditions [10], lo- 
cal spontaneous infiltrations can be achieved. Whether the 
liquid is able to then spread and fill the micro-structures 



depends of the surface patterning. In particular, when 
the surface roughness becomes comparable to an intrin- 
sic value (critical value), a smooth filling is replaced by a 
step-like process in which rows of the structured surfaces 
fill through a zipping mechanism [10]. This gives rise to 
slowing down in the system and the emergence of a dy- 
namics with different time scales. The aim of this paper is 
to characterize this scenario experimentally and with the 
help of numerical simulations. 

For the experiments presented here we choose 
micrometer-sized square pillars set on a regular square 
lattice with height h = 10 /im, width w = 5 /im, spac- 
ing a — 5/im, and hence d = w + a = 10 fim as lattice 
wavelength (see fig. 1(a)). Fabrication of highly precise 
and controllable micro-structured surfaces becomes possi- 
ble through a micro- molding technique [11]. As material, 
we choose a polymer fabricated from a solution of styrene- 
butadiene-styrene (a block-copolymer commercially avail- 
able as "kraton") dissolved in toluene, and then applied 
on a silicon mold prepared by standard photolythography 
and etching techniques. After evaporation, the resulting 
micro-patterned thin polymeric film is bonded to a 170 fj,m 
thick microscope glass slide. Equilibrium contact angle 
for water is of 9 ~ 160 in the metastable CB state and 
8 ~ 142° in the W state (fig. l(b-c)), whereas contact 
angle for a smooth surface made of the same material is 
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Fig. 1: (a) Scanning electron microscopy (SEM) picture of the patterned surface. Here h — 10 /J,m, w = 5 /im, and a = 5 /jm. 
The water droplet spontaneously transits from the metastable CB (b) to the W state (c), resulting in a lower contact angle (AVI 
movie, 0.25 Mb), (d-f) The microstructure underneath the liquid fills in a few hundreds of milliseconds through a square-like 
pattern involving a zipping mechanism (AVI movie, 4.4 Mb) . (g) Enlargement of the zipping mechanism: a row perpendicular 
to the front direction is rapidly filled. The bar in (d-f) indicates 100 fim. 



of 9 f ~ 100°. 

Water droplets are softly deposited on the surface with 
a syringe pump (PHD 2000, Harvard Appartus GmbH, 
March-Hugstetten, Germany) at a very low flow rate 
(5 /il/min) from the outlet of a vertical thin tubing (outer 
diameter 0.158 mm). The outlet is set parallel to the flat 
film, 2 mm above it. A typical droplet grows slowly and 
reaches the dry surface within a minute. Then the flow 
is stopped and a stable drop of contact area ~ 1mm 2 
in the CB state is observed. Observations are performed 
from the bottom of the film with an inverted microscope 
(Axiovert 40 CFL, Carl Zeiss BV, Weesp, The Nether- 
lands). Proper illumination is obtained with a fiber lamp 
(ILP-1, Olympus, Zoeterwoode, The Netherlands), com- 
bined with a high-speed charged coupled device camera 
(APX-RS, Photron Limited, UK). Focusing in the mate- 
rial features is achieved with a piezoelectric objective-lens 
positioning system with reproducibilitybelow 100 nm (MI- 
POS 500, Piezosystem Jena GmbH, Jena, Germany). The 
translucent film is only ~ 40 fim thick, allowing for high 
quality optical imaging of the CB to W transition through 
its bottom. 

The large scale dynamics is captured with a 10 x mag- 
nification objective at 5 kfps. Smaller details are resolved 
at 50 kfps with an oil immersion Plan-Apochromat 100 x 
objective, of numerical aperture NA = 1.4, allowing for a 
few microns thick image plane. 

Submicron resolution is achieved with the help of a 
home-made image post-processing technique. First a 
movie with region of interest 20 /im x 20 y,m is decom- 
posed in an image sequence. For each image, the initially 
117 x 117 pixels window undergoes a lOx oversampling up 
to 1170 x 1170 pixels. Next, a FFT bandpass filter consist- 



ing in a twofold 2d Gaussian filtering in the Fourier space 
is applied. The low-pass component acts as a small struc- 
ture smoother, whereas the high-pass one corrects large 
scale illumination inhomogeneities and increases contour 
contrast. Finally, adaptive intensity threshold is properly 
applied and interface position detected. Attention should 
be paid on the fact that this method is sensitive to illu- 
mination variations through the image sequence, and may 
result from time to time in a slight shift of the detected 
interface position. Moreover, shadowing caused by pillars 
yields an imprecise detection of the interface position in 
the vicinity of the pillars, although the main part is accu- 
rately detected within 100 nm. 

After depositing the drop on the surface, a sponta- 
neous breakdown of the metastable CB state starts usually 
within a minute. From a local infiltration point, a rapid 
spreading (responsible for the lateral front, hereafter called 
"main front") develops through the square lattice posts. 
The details of this local infiltration point are sensitively 
dependent on the height (i.e. roughness): the higher the 
posts, the larger the time needed to achieve that. This 
is probably due to the presence of some local energy bar- 
rier [9] whose value is increasing with surface roughness. 
In this way, for large h, the infiltration threshold cannot 
be achieved spontaneously and some local trigger mecha- 
nism is therefore necessary. Once the liquid has entered 
the micro-structures, the dynamics of liquid moving within 
the pillars can be understood in terms of energy balance 
between the pulling mechanism of the liquid-gas interface 
on the top of the pillars (energy gain) and the repulsive 
wetting of the hydrophobic walls (energy loss). Therefore, 
as shown in [10], whether the front propagates depends on 
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v z (mm/s) Side length (jam) 

Fig. 2: (a) Temporal evolution of the front position. Values have been averaged over the 4 sides of the square-shaped pattern, 
starting from its center, (b) Probability distribution of the zipping velocity V z . (c) Measured latency as a function of the side 
length. How latencies have been calculated is explained in the text. Acording to the theory, the latency becomes negative at 
a side length of 270 ^tm (see arrow). As observed in the corresponding movie, beyond this value the pattern is not a perfect 
parallelepiped anymore. 



the ratio a/h. More precisely, if < C with 



then filling can occur. The contact angle Of for water on 
flat kraton is Of « 100°, thus very close to the critical con- 
tact angle C = 101.5° set by a — 5 fim and h = 10 fj,m. 
Therefore, the dynamics is slow and a square-shaped pat- 
tern grows in a stepwise manner, by successively filling 
the rows along the 4 sides (zipping regime, fig. l(d-g)). It 
has also been pointed out in [10] that this critical slowing 
down (involved when Of becomes close to C ) results in 
large dispersion of the measured front velocities, carried 
out from several experiments with the same sample, al- 
though in each case the front velocity remains almost con- 
stant during the filling. It can be noticed that, although 
averaged front velocities vary from one experiment to the 
other close to the critical point, the main features persist, 
provided that experiments are performed close enough to 
the critical point (i.e. the zipping persists). The ques- 
tion we address in this letter is: what characterizes the 
dynamics during the invasion process? 

To answer this question, we have first performed large 
scale measurements. Figure l(d-f) shows snapshots dur- 
ing filling while fig. 2(a) shows the main front position 
as a function of time in a typical experiment. When the 
filled area only extends over few pillars the dynamics un- 
dergoes a transitional regime, also observed in the numer- 
ical simulations. Afterwards, the front velocity is found 
to be Vf = 2.3 mm/s. Since Vj is almost constant and, 
at the same time, the square-shaped pattern remains, one 
may wonder how the zipping process adjusts. Figure 2(b) 



shows the probability distribution of the zipping velocity 
V z involved to fill a row, throughout the filling, processed 
over approximately 500 frames. Interestingly, it shows a 
bimodal distribution with a significant dispersion. One 
may speculate that the maximum of the probability dis- 
tribution, located at y™ ax = 62 mm/s, corresponds to the 
wetting dynamics for an ideal sample. The large disper- 
sion can be attributed to surface inhomogeneities resulting 
in a surface energy lowering (caused by defects or impu- 
rities). Larger defects might be the cause of the second 
peak, located around V z — 105 mm/s. 

Beyond the observed dispersion of the zipping velocity, 
it is worthwhile noticing that rows do not fill one by one 
continuously. A schematic picture of this is provided in 
the inset at fig. 2(c). When the filled square is still small, 
a sideway row gets filled (at time T±), and then the next 
row (at time T 2 ) after a latency time (It = T% — T\ > 0). 
The latency is decreasing with increasing square size until 
it finally vanishes (at about 270 [im square side length, ® 
in fig. 2(c)). Afterwards, a different scenario is observed. 
The new row is already started (at time T{) before the for- 
mer one has been completed (at time T' 2 \ In this case, we 
express the latency time as k =T[—T' 2 < (© in fig. 2(c)). 
As a consequence, the filled area cannot be seen as a per- 
fect parallelepiped any longer although the macroscopic 
square-like shape remains, considering that the filled sur- 
face gets larger (see fig. l(d-f)). Moreover, the side length 
corresponding to l t — can simply be calculated from 
velocities by V™ ax /Vf x d = 270 fim, which is in good 
agreement with the length extrapolated in fig. 2(c). 

Details of the front advancing in between the posts are 
shown in fig. 3(a). The height is a priori not exactly 
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Fig. 3: (a) Details of the zipping process in the experiment. Liquid-gas interface positions obtained every 20 fis are simultaneously 
drawn. The filling begins at bottom-left corner. The main front (responsible for the macroscopic square shape) is from bottom 
to top. The zipping front (in a row to be filled) is from left to right. Color changes after an interface has detached from a pillar; 
in chronological order: blue, red, green, pink, orange (AVI movie, 2.6 Mb), (b) Interface position as a function of time along 
the 4 straight lines suggested by the arrows in (a). Curve and arrow colors match. Labelled rectangles correspond to different 
steps in the dynamics, as explained in the text. 



known. The dynamics of the interface are captured in a 
plane close to the middle of the pillars. Based on the ob- 
jective specifications, we estimate it to be at z = 4 ± 1 \im 
from the base. The dynamics are processed from an ex- 
perimental record chosen especially for its low V z /Vf ratio 
to better reveal main and zipping front progressions at 
once. Thus, when averaged over 20 /im, the main front 
velocity is Vf = 1.9 mm/ s and the zipping front velocity 
is V z — 10 mm I s. The figure enables a comprehensive pic- 
ture of the complex dynamics in the 2d plane: position, 



curvature and velocity variation of the interface are easily 
accessible. Front and zipping interface positions as a func- 
tion of time are quantitatively tracked in fig. 3(b) along 
the straight lines passing through the middle of the 4 rows 
shown in fig. 3(a). First, the two perpendicular liquid- 
gas interfaces defined as main and zipping fronts spread 
between the pillars and reach almost simultaneously the 
edge. Next, they fastly merge, and the resulting inter- 
face straighten. Rectangles drawn in fig. 3(b) highlight 
this stage. In rectangles labelled "1", main and zipping 
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Fig. 4: Numerical simulations reveal the details of the zipping mechanism. Two values of the wetting properties are used with 
respect to the critical angle given in equation (1): cosd/ cos# c = 0.868 ((a) and (b)) and cos#/cos# c = 0.927 ((c) and (d)). 
The interface position is tracked as a function of time (from left to right) at a given height z (z/h = 0.3 for (a) and (c) and 
z/h = 0.3 for (b) and (d)) and all plots are equispaced with the same numerical time lag. When approaching the critical point 
((c) and (d), (AVI movies, 0.2 Mb)) the zipping is experiencing a dynamical slowing down followed by a rapid merging process. 
The pulling mechanism provided by the liquid on the top of the pillars is inducing a complex 3d stretching that deforms the 
interface. An insight of this is provided by the cuts at different heights (z/h = 0.5 and z/h = 0.3). In particular, quantitative 
agreement is found with the experiments, assuming a height less than h/2 (z/h = 0.3 in (c)). 



fronts are superimposed, indicating that the interface at 
an intersection is almost symmetric. In rectangles labelled 
"2" , the merged front has reached a lower velocity and an 
asymmetry develops: V z increases more than Vf. After- 
wards, when the interface has reached the opposite post, it 
splits into 2 fronts that progress in main and zipping direc- 
tion respectively, with drastically different velocities. The 
passage is made within 4 ms for the main front whereas it 
is almost 2 orders of magnitude less for the zipping front. 

Here the interface is described in one plane, but since 
the filling is a 3d process, it is straightforward that sev- 
eral planes should be analyzed simultaneously. Unfortu- 
nately this cannot be achieved experimentally. Therefore, 
we have carried out numerical simulations to better reveal 
the 3d structure as well as the dynamics of the zipping 
mechanism. 

For the numerical simulations presented here, we use 
a 3d lattice Boltzmann [12] model for single component 
multiphase flows whose details are discussed in [13]. In 
all these mesoscale models [14], it is impossible to match 
all the physical parameters and preserve at the same time 
a feasible computation. Typically the width of the inter- 
face is too large and the liquid-gas density difference is 
unphysically small with respect to physical reality. As al- 
ready noticed elsewhere [15], this may affect the speed of 
the interface that, if directly proportional to the interface 
width and inversely proportional to the density difference, 
will move too quickly in the simulations when compared 
to an experiment. This must be taken into account by 
renormalizing the time scale if one wants to reconstruct a 
correct pathway as a function of time [15]. 

All the geometrical aspect ratios for the surface micro- 
structures can be adapted in order to match those typi- 
cally studied in the experiments [16] and wall interactions 
with wetting properties are introduced as explained in [13]. 
The initial condition consists of a liquid reservoir on the 



top of the pillars and an infiltration characterized by a flat 
front perturbed with an advanced liquid precursor, just as 
experimentally observed at the onset of zipping. In this 
way the liquid precursor induces a lateral motion that is 
continuously driven by the pulling mechanism of surface 
tension on the top of the pillars. After a small time tran- 
sient, the properties of this numerical zipping motion show 
a remarkable independence on the initial condition details 
and are reported in fig. 4 (zipping is advancing from left to 
right in these figures). Two wetting properties are shown, 
cos0/cos0 c = 0.868 ((a) and (b)) and cos6>/cosfl c = 0.927 
((c) and (d)) and the interface contour is tracked in a two 
dimensional plane for two different heights, z/h = 0.3 ((a) 
and (c)) and z/h = 0.5 ((b) and (d)). 

The results presented in figs. 3(a) and 4 reveal the com- 
plex structure of the zipping process, both from experi- 
mental and numerical viewpoints. Firstly the liquid pre- 
cursor fastly merges with its neighbor and after that expe- 
riences a passage through the walls. It is observed that as 
one approaches the critical point the numerical pathway 
at the height z/h — 0.3 in fig. 4(c) shares quantitative 
agreement with the experimental one. The main steps of 
the zipping process are sketched in fig. 5 and should be 
superimposed to the above figures. 

Starting with the case z/h = 0.3, the successive steps 
are described as numbered in fig. 5. 

1- The pulling mechanism induced at the top of the 
pillars makes zipping and main front interfaces merge (red 
point at the corner of pillar A). 

2- The liquid bridge suspended on the four pillars col- 
lapses at pillar A, producing a local top-to-bottom pulling 
mechanism. The dynamics is dominated by surface ten- 
sion effects that tend to reduce the strongly curved surface 
of the merged interface, by drawing liquid from the top. 
Wall effects are negligible in this extremely fast process 
(within 20 us in the experiment): it is almost independent 
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Fig. 5: Sketch of the different steps in the zipping process. The 
numbered areas refer to the text. As in figs. 3(a) and 4, the 
zipping front is from left to right (green arrows) and main front 
is from bottom to top (red arrows). 

of the wetting properties (compare for example case (b) 
and (d) in fig. 4). 

3- Since the interface is straight, the pulling persists but 
the interface above is behind (late). Viscous dissipation 
dominates, making the front much slower (time scale is 
« 1 ms in the experiment). Moreover, the interface is still 
pinned on the edges of pillars B and C. 

4- The interface reaches pillar D and splits into two parts 
and then the filling progresses. The dynamical details are 
influenced by wall interactions that cause a slowing down 
for the propagation. This is strongly pronounced espe- 
cially when the wetting properties approach the critical 
point passing from (a) to (c) in fig. 4. 

The pulling mechanism on the top of the pillars is 
stretching the interface, and obviously the closer to the 
top the stronger it is. Consequently, we choose a cut plan 
at z/h = 0.5 to evidence the change in the dynamics. If 
we compare with z/h — 0.3, the front is slower in region 
2, but much faster in region 3 where the slowing down is 
not observed (compare (c) to (d) in fig. 4). Interestingly, 
this makes the liquid invade the gap between pillars C and 
D at z/h — 0.5 before than at z/h = 0.3. The resulting 
stretching overcomes the wall interactions that, being in 
the slightly hydrophobic regime (Of ~ 100°), act to pre- 
vent a complete wetting of the boundaries (figs. 3a and 
4c). 

One has to remember that liquid is present in area 5, 
causing the interface that bridges over pillars C and D to 
be more curved than the one over pillars B and D (the 
portion beyond these pillars in not filled). This induces 
a faster filling in the zipping direction (green arrows in 
fig. 5) than in the main front direction (red arrows). The 
front velocity between pillars C and D is progressively in- 
creased while the interface along pillar C is deformed by a 
stronger pulling (fig. 4(d)). It can be assumed that this ef- 
fect is mainly at the base of the macroscopic square front 
and that the zipping velocity is quite sensitive to slight 
variations of the shape of the pillars top. 

In summary, we have investigated experimentally and 
numerically details of the spreading dynamics for the CB 
to W transition on a micro-patterned substrate, close to 
the critical point above which filling occurs. The already 
wet portion expands by filling the rows along the pattern 



through a zipping mechanism, thus preserving at a macro- 
scopic scale the structure of the underlying square lattice. 
We have shown that in this regime, even if the pattern 
grows in a regular way (the front velocity Vt is almost 
contant), the zipping velocity V z involved to fill a row 
undergoes a noticeable dispersion (from 25 to 125 mm/s 
in a typical experiment). This reveals the strong sensi- 
tivity of the dynamics to local wetting properties with- 
out prejudging the cause (roughness variation, geometrical 
imperfection, material inhomogeneity, dust deposit, ...). 
Interestingly, this translates into a latency acting as a self- 
adjusting parameter and, consequently, makes the square 
shape robust at a macroscopic level. Interested in getting 
further details at a smaller scale, we have performed, to 
our knowledge, unprecedented high-speed submicron ob- 
servations of the zipping- wetting dynamics. The dynami- 
cal details of the interface demonstrate the complexity of 
the invasion process, while being in very good agreement 
with numerical predictions for z/h — 0.3 (corresponding 
to z = 3/zm in the experiment). Moreover, the numerical 
results shown for z/h = 0.5 give a deeper insight of the 
complex 3d stretching experienced by the interface. Fi- 
nally, the analysis of the submicron dynamics allows us to 
provide a simplified explaination for the resulting wetting 
dynamics. 
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